Zero-energy states of massive Dirac equation in magnetic fields 
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The Dirac equation with a U{1) vortex in the mass-term is solved in the presence of magnetic-like 
fields at zero energy. By drawing an analogy to classical mechanics it is shown that the four- 
component Dirac equation in arbitrary magnetic field always yields one zero-energy state. In the 
time-reversal preserving, pseudo-magnetic field, however, the number of zero-energy states may 
depend on the field's profile and sign. Some explicit examples are worked out. Possible implications 
of these results for the electrical charge of the vortex and for the behavior of graphene in magnetic 
field are discussed. 



I. INTRODUCTION 

Relativistic invariance dictates that the spectrum of 
Dirac's equation extends over both positive and negative 
values. Even in the condensed-matter context where the 
invariance is only pseudo-relativistic, the spectrum is of- 
ten symmetric with respect to the change of the sign of 
the energy. Particularly interesting then are the states 
with precisely zero energy, right in the middle of such a 
spectrum. They have been known to exist in at least two 
sets of circumstances: when the Dirac electron is sub- 
jected to the magnetic field or when the mass-term 
forms a kink or a vortex in the configurational space Q . 
In either case the number of the zero-energy states is 
related to some global property of the external potential 
that acts on the electron, such as the magnetic flux or the 
vorticity, but is otherwise independent of the details of 
the configuration. The existence of the macroscopically 
degenerate manifold of zero-energy states is believed to 
be responsible for some of the observed quantum Hall ef- 
fects in graphene in uniform magnetic field Q , for exam- 
ple. Vortex zero-energy states have been intensely stud- 
ied theoretically for a number of years, within a scenario 
for the fractionalization of electric charge Q • In the con- 
text of graphene, in particular, the zero-energy states in 
the core of the vortex of a superconducting order param- 
eter have recently been argued to provide an example of 
the elusive Majorana fermion [H, as well as to render the 
core of the vortex ordered @. 

It is less clear if, and what kind of zero-energy states 
exist in presence of both the vortex in the mass-term in 
the Dirac equation, and the magnetic fields. This prob- 
lem arises in several current lines of investigation. In 
the context of fractionalization, Jackiw and Pi [7| have 
shown that the addition of a localized flux of the pseudo, 
that is time-reversal symmetry preserving, magnetic field 
does not change the number of zero-energy states, but 
only modifies their form. The problem of the vortex in 
the superconducting order in graphene can be mapped 
onto the same Hamiltonian [6]. The recent observation 
of the Kosterlitz-Thouless scaling of the resistivity near 
the metal-insulator transition in graphene also urgently 
calls for a better understanding of the vortex structure of 
Dirac fermions in true, time-reversal symmetry breaking, 
magnetic field [8]. At the present time there are at least 



two candidates for the possible order parameters in this 
system with the requisite U{1) symmetry: the Neel or- 
der parameter, favored by the on-site repulsion between 
electrons, which acquires an easy plane in the magnetic 
field ti] , and the Kekule bond-density- wave (BDW) order 
parameter, which may be favored by the electron-phonon 
coupling fl3, im . The internal structure of the vortex in 
either of the two order parameters in vanishing magnetic 
field has been a subject of several investigations in the 
recent past [lol - [ll |. 

Here I consider the problem of zero-energy states of the 
two-dimensional Dirac Hamiltonian with a unit vortex in 
the mass-term, in the presence of either pseudo or true 
magnetic field. There are two principal results. For the 
pseudo-magnetic field, a sufficiently extended field's con- 
figuration may change the number of zero-energy states 
to zero or two, depending on the field's direction. In case 
of a finite pseudo-magnetic flux, on the other hand, the 
number of zero-energy states remains at unity, in agree- 
ment with previous work Q. For true magnetic fleld, 
the Dirac equation at zero-energy is shown to be mathe- 
matically equivalent to classical Hamilton's equations in 
a certain time-dependent potential. This analogy is used 
to show that there is always precisely one zero-energy 
state irrespectively of the field's profile. This agrees with 
the solution of the special case of a uniform mass and 
a field The resulting zero-energy state is naturally 
time-reversal asymmetric, with the asymmetry growing 
with the magnetic field. Finally, it is argued that the vor- 
tex in a general U(l) order parameter in graphene may 
carry the electrical charge of zero or one, depending on 
the type of orders both outside and inside the vortex core. 
I close with the short discussion of the recent experiment 
in light of these results. 



II. EQUATION FOR ZERO-ENERGY STATES 

We are interested in the spectrum of the general two- 
dimensional Dirac Hamiltonian in presence of the vortex 
configuration in the two-component mass and the vector 
potential: 

H = ijojiiPi - CAi{f)) - Ai(f)?;7o73 - A2(f)i7o75, (1) 
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where i = 1,2, the summation over the repeated indices is 
assumed, and the matrix C G {/,735}, with 735 = 17375. 
We chose the vortex configuration to be rotationally sym- 
metric, sothat Ai(f) = |A(r)|cos6', A2(r) = |A(r)|sin6', 
where {r, 9} are the polar coordinates. The 7-matrices, 
as usual, form a representation of the Clifford algebra, 
and [7i,7j]+ = 2^^, for i = 0,1,2,3,5. 

Although the results will be representation indepen- 
dent, to be specific we will use the representation directly 
relevant to graphene [l^, in which ^70 71 = — 173 (81 Ci, 
«7o72 = -/2«)cr2, ■i7o73 = -(Ji<E)(Ji, and 17075 = cr2(S)<Ji. 
{/2,(t} is the Pauli basis. The first Pauli matrix acts in 
the valley, and the second in the sublattice space. In 
this representation 70 = I2 ® (^Zi and it anticommutes 
with the Hamiltonian. In graphene, the Hamiltonian 
in this representation would act on the Dirac fermion 
^I' — M_, with u± {v±) representing the 

wavefunction on sublattice A (B) with Fourier compo- 
nents in the vicinity of the Dirac point ±K. When C = I 
the vector potential represents the true, time-reversal 
symmetry breaking magnetic field, and when C = 735, it 
stands for the time-reversal preserving pseudo-magnetic 
field, which would represent the low-energy effect of rip- 
ples of the graphene sheet, for example [1^]. Note that 
in the above representation 735 = (J3 /2 , so that in the 
latter case the sign of the coupling of the magnetic field 
to Dirac fermion is different for the two Dirac points. 

The specific mass-term in the graphene representation 
presented above would describe the vortex configuration 
in the BDW. The spectrum is nevertheless completely 
general, as any Hamiltonian with the form of Eq. (1) 
can be transformed into the graphene representation by 
a unitary transformation [6]. This follows from the fact 
that the four anticommuting matrices appearing in the 
Hamiltonian form a representation of the Clifford alge- 
bra, and all such representations are equivalent [17]. 

Since the zero-energy states are the eigenstates of 70, 
it is useful to transform to the representation in which 
7o would become the block-diagonal matrix a^® l2- In 
the graphene representation this is easily accomplished 
by exchanging the two sets of Pauli matrices. The 
Hamiltonian then becomes purely off-diagonal: H ^• 
(tT+ ® Hy + a- ® i?„)/2, with (t± = cti ± ztT2, and the 
equations for the zero-energy states are: 

i{2d-, + ^z)u+(f) + A(f)u_(f) = 0, (2) 
r 

- i(2d, - s^z)u-(r) + A(f)u+(f) = 0, (3) 
r 

with v-f- = V- = 0. Here, z = x + iy, dx + idy — 2dz, 
and A = Ai + iA2. We are using the symmetric gauge 
in which A{r} — x(^)(~yj 2;), so that the field strength is 
rotationally symmetric and determined by B(r) = (dr + 
{l/r))x{r)- The sign s = 1 refers to the real, and s = — 1 
to the pseudo-magnetic field. 

Note that the other possibility, namely u± = 0, H^v — 
0, does not yield a normalizable zero-energy state. This is 



because assuming a finite, real or pseudo, magnetic field 
at the origin, ''j ^-^id it can be neglected. This 

implies that v{r) cx 1/r near the origin, just as without 
the magnetic field, and one must choose v±(r) = 0. For 
the antivortex the role of u and v are reversed. 

Rotating the variables and assuming rotational sym- 
metry of the solution, as u+(r) = \/^^q{r) and u_(r) — 
Vip{r), we may finally write the equations as 

drq{r) = -\AirMr) - x{r)q{r), (4) 

drPir) = -\A{r)\q{r) + sx(r)p(r), (5) 
which are solved in the remainder of the paper. 

III. PSEUDO-MAGNETIC FIELD 

Consider the simpler case of the pseudo-magnetic field 
with s — —1 first. By rotating in the (g,p)-plane by 7r/4 
it immediately follows that the general solution is 

p{r)±q{r) - ^^g- /;(±|AK)l+x(r-'))*r'_ (g) 

Assuming that the amplitude |A(r — > 00 )| — > Aq, when 
x{r) = normalizability of the state dictates that the 
constant A- — 0, and therefore p{r)_= q{r). In this case 
there is a single zero-energy state 2]. If xi^) 7^ 0, but 
\x{r — >■ 00) I < Ao, it is still A- — 0, and there is still a 
single zero-energy state. If xi^) < Oj but \x{r — > 00) | > 
Aq, there are no normalizable zero-energy states at all, 
as both constants A± — 0. 

When x(^) > Oj ^-i^d ^ 00) ^ with a > 0, 
on the other hand, there are evidently two orthogonal 
normalizable zero-energy states. We may write them as 

4=A/-±e-/o'^(±''^('-')l+^('-'»'^'^'(±l,z), (7) 

where Af± are the normalization constants. The norm of 
both states will be finite if the strength of the field B{r) 
decays slower with radius than 1/?'. In particular, if the 
pseudo-magnetic field and the amplitude of the mass are 
uniform, whereas the first state is centered at the origin 
the second state is sharply peaked at r^ax — 2Ao/B, 
which diverges as the field approaches zero. 

In sum, for a sufficiently localized pseudo-magnetic 
field of either sign, there is precisely one zero-energy state 
localized at the origin Q. If the field decreases slower 
than the inverse radius, on the other hand, the number 
of zero-energy states is zero or two, depending on the 
field's direction. 



IV. MAGNETIC FIELD: ANALOGY WITH 
MECHANICS 

For the true magnetic field (s = 1), assuming |A(r)| 
to be monotonically increasing, we may introduce a time- 
like dimensionless variable t as dt = \A(r)\dr, and rewrite 
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Eqs. (4)-(5) in a more suggestive form as 

Pit) = -q{t) + f{t)p{t), 



(8) 
(9) 



where x = dx/dt, and f{t) — In this form 

the Dirac equation for the zero-energy state may be rec- 
ognized as Hamihon's equations of classical mechanics 
for the "coordinate" q and the "canonical momentum" 
p, with the time-dependent classical Hamilton's function 



Hiq,p,t) = -f{t)pq+^{q^-p^) 
and with the corresponding Lagrangian 

^L{q,q,t) - ^q^+'^q\ 



(10) 



(11) 



with K{t) = 1 + f{t)^ — f{t), and with a total time deriva- 
tive omitted. Apart from the irrelevant minus sign, this 
Lagrangian describes the motion of a classical particle in 
the inverted harmonic potential 



Viq,t)^-^-^q^. 



(12) 



In absence of the magnetic field K{t) = 1, and the 
potential V{q,t) is static; for any initial q{0) providing 
the right amount of initial kinetic energy will get the 
particle to the origin in infinite time. Since the initial 
energy is simply E(t = 0) = (p^(0) — q'^{0))/2, this is 
ensured by the initial condition p{0) ~ g(0). The con- 
servation of energy implies then that p{t) — q{t) at all 
times. This mechanical interpretation of the familiar so- 
lution becomes useful in understanding the qualitative 
effect of the magnetic field. Although for a finite mag- 
netic field the potential V{q, t) typically becomes steeper 
with time, it is clear that there are two independent so- 
lutions of the equation of motion: first when the initial 
kinetic energy is too small so that particle starting from 
some g(0) > returns and ultimately runs off to posi- 
tive infinity, and the second, when for a too large initial 
kinetic energy the particle goes over the top, and runs 
off to the negative infinity. Continuity guarantees then 
the existence of the initial condition in between these two 
extremes for which the particle will reach the top with 
precisely zero velocity, i. e. in infinite time. This trajec- 
tory corresponds to the normalizable zero-energy state 
of the Dirac Hamiltonian. Using the equation of motion 
the rate of the change in the " mechanical energy" of the 
system is found to be 



E = -^qHt). 



(13) 



To reach the top with zero velocity the sum of the ini- 
tial energy and the work done externally needs to vanish. 
This may be expressed as a global condition on the solu- 
tion. 



n{t)q{t)q{t)dt = 0. 



(14) 



V. EXAMPLE: UNIFORM MAGNETIC FIELD 

For illustration, let us consider the case of the uniform 
magnetic field x(^) = Br/2, and the vortex with the core 
of size R: |A(r)| = Aor/i?, for r < R, and |A(r)| = Ao, 
for r > R. The equation of motion is 



{l + f{t)-f{t))q. 



(15) 



For r < R, f{t) — BR/2Aq and constant, and t 
Aqt^ /2R. The general solution is therefore 



g<(r) cie^ + C2e 



(16) 



where Ci and C2 are constants, and L is the characteristic 
length, determined by 



L- 



(17) 



When r > _R, on the other hand, / = St, where S — 
B/{2{/S.qY) is a dimensionless parameter. It would then 
seem that one needs to distinguish two cases: 
a) when 5 < 1, 



q>{r) = C3(Aor)e" 



■u 



. i\5\ + \-5 3 \By 



(18) 



where U stands for the hypergeometric function. In this 
case the potential V{q, t) is always repulsive and conse- 
quently the solution is monotonically decreasing, 
b) when 6 > 1, the solution is 



'7>(^) = C3(Aor)e~ 



45 



25 3 \B\r^ 
2' 2 ^ 



(19) 



Although the potential V{q,t) now starts as attractive, 
the solution is nevertheless still always monotonically de- 
creasing, and without oscillations. It is in fact qualita- 
tively the same as in the previous case. 

The continuity implies that g<(i?) — q>{R) and 
dq^{R)/dr = dqy{R) / dr, which yields two linear equa- 
tions on the three constants q, i = 1,2,3. The normal- 
ization then provides the third condition that completes 
the solution, as usual. It is also easy to see that the 
solution is always with C2 > ci and therefore monoton- 
ically decreasing, as one indeed would expect from the 
mechanical analogy. 



VI. REMARKS ON THE SOLUTION 

The reversal of the sign of the true magnetic field sim- 
ply exchanges the solutions for p and q, as evident from 
Eqs. (8)-(9). This is of course equivalent to the time re- 
versal, consisting of the exchange of the two Dirac points 
followed by the complex conjugation [l5| . 

The solution for the antivortex may be obtained most 
simply by multiplying the zero-energy state by the ma- 
trix 75, for example. Since 75 anticommutes with the 



4 



last term in Eq. (1) while commuting with the rest, the 
Hamiltonian ^^H^^ has only the sign of A2(r) flipped, 
and thus represents an antivortex. In the graphene rep- 
resentation 75 = (72 (81 172, and the result of the multipli- 
cation correctly reproduces the time reversal followed by 
the exchange of the sublattices. 

It is interesting to consider the two limits of the above 
zero-energy solution, for weak and strong magnetic field. 
In the former case the zero-energy state approaches the 
standard zero-field solution of Jackiw and Rossi which 
respects the time-reversal symmetry. In the latter case, 
when B ^ Ao/i?, the solution for q{r) is localized very 
close to the center of the core, and thus it has the form 
of the second term in Eq. (16). Eq. (4) then fixes the 
second component to be p{r) ~ {BR/AQ)q{r) ^ q{r). In 
the limit of strong magnetic field the zero-energy solution 
is, in the graphene representation. 



0) + 0(Ao/Bi?), (20) 



and has a single large component. Note how the result 
does not approach continuously the solution at A = 0, 
since even for an infinitesimal A vorticity eliminates the 
w-components of the state completely. When the mag- 
netic field is much larger than the mass the zero-energy 
state is strongly localized, finite only over one sublattice, 
and has Fourier components essentially only near one of 
the Dirac points. 

Since the Hamiltonian in Eq. (1) always anticom- 
mutes with 70, even a rotationally asymmetric mass- 
vortex and/or the magnetic field still produce a spectrum 
symmetric around zero. This means that continuous de- 
formations of the rotationally invariant problem consid- 
ered here would still yield exactly one zero-energy state, 
because moving that single state to a finite energy would 
obviously violate the spectrum's symmetry. 



VII. ELECTRICAL CHARGE OF THE VORTEX 

The existence of the zero-energy solutions in the mag- 
netic field being established, let us turn to some of their 
possible manifestations. Let us assume first that the or- 
dered state of the electrons in the magnetic field at half 
filling, i. e. at the filling factor zero, preferred by the 
electron-electron interactions is indeed the BDW [Tcil.[Tl|. 
identical for both projections of the spin, for example. 
This being an order parameter with an (approximate) 
U(l) symmetry, a vortex will be topologically distinct 
and thus stable. The mean-field Hamiltonian in presence 
of such a vortex will then consist of two (spin) copies of 
the Eq. (1) in the true magnetic field, which will have two 
zero-energy states. The zero-energy Hilbert space then 
provides a representation of the algebra of the possible 
order parameters in the core Q . In the case at hand this 
algebra consists of the arbitrary oriented Neel and the 
charge-density-wave (CDW) order parameters, given by 
{cj (g) 7o,/2 ^ 7o}, respectively, since these four matrices 



exhaust the set of operators that anticommute with the 
Hamiltonian I2 ^ H- The electrical charge of the vortex 
vanishes if the Neel order parameter is preferred locally, 
since that implies that precisely one state from the zero- 
energy space is occupied and the other orthogonal state 
is empty. If, in contrast, the CDW state is preferred en- 
ergetically, both states will be occupied (for vortex) or 
empty (for antivortex) and the vortex will bind a unit 
charge [iJl . Obviously there is a cost in the electrostatic 
energy in the latter case, but one can conceive of, possibly 
unrealistic, Hamiltonians, with the nearest-neighbor suf- 
ficiently strong, for example [l^, where such a state could 
be energetically favorable. At finite temperatures, fluctu- 
ating vortices and antivortices would in that case indeed 
behave like a plasma of positive and negative charges, 
with the metallic conductivity proportional to the vortex 
density. Assuming that density vanishing at the point 
of metal-insulator transition in the standard Kosterlitz- 
Thouless fashion would then rationalize the experimen- 
tal observation P. This scenario, however, obviously 
depends on the favorable ordering of the interaction en- 
ergy scales in the many-electron Hamiltonian, so that 
electron-phonon coupling is larger than nearest-neighbor 
repulsion, which in turn is larger than the on-site repul- 
sion, for example. Recent work [l^ which finds the stable 
BDW over a large portion of the zero-field phase diagram, 
however, could make this mechanism more plausible. 

The second option mentioned in the introduction for 
the order parameter with the (exact) U(l) symmetry is 
the familiar Neel order, which for small enough Zeeman 
energy develops an easy-plane @. Vortex in this order 
would be described by the mass-term in the Eq. (1) re- 
placed by (AifTi -I- A2cr2) ^ 7o- This being equivalent to 
Eq. (1) , the Hamiltonian for the Neel vortex also contains 

This time 
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two zero-energy states m its spectrum 
however, the algebra of competing orders in the vortex 
core is different, and is readily seen to consist of the third 
component of the Neel, two BDWs, and the third compo- 
nent of the Haldane-Kane-Melle (HKM) order parameter 
m, {^s 70, h ® ilol?.,h ® ilolb, 0-3 17172}, respec- 
tively. In this case the vortex will be charged only if it is 
the HKM ordering that is energetically preferred in the 
core, and otherwise not. In particular, one would assume 
that the order in the core, if possible, should be of the 
same type as the order supporting the vortex, which was 
preferred energetically in the first place. On this basis it 
is the third component of the Neel order that may be ex- 
pected to develop in the core [l^, with the consequence 
of rendering the vortex electrically neutral. 



VIII. CONCLUSION 

In this paper it has been established that the two- 
dimensional Dirac equation with a vortex in the two- 
component mass-term has a single zero-energy solution 
even in presence of an arbitrary, including infinite, flux 
of the (true) magnetic fleld. When the particle in ques- 
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tion has spin-1/2, the number of zero-modes is doubled. 
In the context of graphene this should imply that, de- 
pending on the residual terms in the Hamiltonian such 
as electron-electron interactions, non-linearities of the 
electron dispersion, edges, etc., the electrical charge of 
the vortex excitation in some C/(l) order parameter in 
magnetic field may be zero or one, corresponding to the 
occupancy of both, or just one of the two zero-energy 
modes. Two candidate order parameters and possible 
connections to the experiment on metal-insulator transi- 



tion in graphene were discussed. 
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